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CONFIRMATION OF HOSTLESS TYPE IA SUPERNOVAE USING HUBBLE SPACE TELESCOPE IMAGING 

M. L. Graham 1 , D. J. Sand 2 , D. Zaritsky 3 , C. J. Pritchet 4 

ABSTRACT 

We present deep Hubble Space Telescope imaging at the locations of four, potentially hostless, 
long-faded Type la supernovae (SNe la) in low-redshift, rich galaxy clusters that were identified in the 
Multi-Epoch Nearby Cluster Survey. Assuming a steep faint-end slope for the galaxy cluster luminosity 
function (ad = —1.5), our data includes all but < 0.2% percent of the stellar mass in cluster galaxies 
(< 0.005% with ad = —1-0), a factor of 10 better than our ground-based imaging. Two of the four 
SNe la still have no possible host galaxy associated with them (Mr > —9.2), confirming that their 
progenitors belong to the intracluster stellar population. The third SNe la appears near a faint disk 
galaxy (My = —12.2) which has a relatively high probability of being a chance alignment. A faint, 
red, point source coincident with the fourth SNIa’s explosion position (My = —8.4) may be either 
a globular cluster (GC) or faint dwarf galaxy. We estimate the local surface densities of GCs and 
dwarfs to show that a GC is more likely, due to the proximity of an elliptical galaxy, but neither can 
be ruled out. This faint host implies that the SNIa rate in dwarfs or GCs may be enhanced, but 
remains within previous observational constraints. We demonstrate that our results do not preclude 
the use of SNe la as bright tracers of intracluster light at higher redsliifts, but that it will be necessary 
to first refine the constraints on their rate in dwarfs and GCs with deep imaging for a larger sample 
of low-redshift, apparently hostless SNe la. 

Subject headings: supernovae, galaxies: clusters 

1. INTRODUCTION 

Supernovae of Type la (SNe la) are the thermonu¬ 
clear explosions of carbon-oxygen white dwarf (CO WD) 
stars, commonly used as cosmological standard candles 
although their progenitor scenario is not yet well under¬ 
stood (e.g., Howell 2011). Most likely, the WD is in a 
binary system with either another WD or a red giant 
or main sequence star, and the explosion occurs after 
merger with, or sufficient mass accretion from, the com¬ 
panion. The explosion rate of SNe la is correlated with 
galaxy mass and star formation rate, and most of the dis¬ 
covered SNe la reside in large galaxies - but since SNe la 
are very bright, they are also used as “cosmic light¬ 
houses” for faint or diffuse astrophysical structures that 
are difficult to assess directly. The purpose of this work is 
threefold: (1) to address the utility of SNe la as a bright 
tracer of baryons in rich galaxy clusters, (2) to confirm 
that SNe la progenitors include truly old (> 2 Gyr) pro¬ 
genitor systems, and (3) to investigate constraints on the 
SN la rate in faint hosts such as dwarf galaxies and glob¬ 
ular clusters (GC). We motivate each of these science 
goals in turn. 

1.1. Baryon Accounting 

Understanding the growth of structure in the universe 
requires a full accounting of baryonic mass, and this must 
include the population of intracluster (IC) stars that were 
stripped from their host galaxy and reside in the clus¬ 
ter potential. Direct measurements of this low surface 
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brightness component are possible (e.g., Gonzalez et al. 
2005; Zibetti et al. 2005; Montes & Trujillo 2014; De- 
Maio et al. 2015), but are difficult beyond the local uni¬ 
verse due to cosmological surface brightness dimming. 
Indirect measurements of the fraction of intracluster light 
(/icl) can be made using bright tracers of the under¬ 
lying stellar population such as planetary nebulae and 
novae, which has been done for the nearby Virgo and 
Fornax clusters respectively (Feldmeier et al. 2004; Neill 
et al. 2005). At higher redshifts a brighter tracer is re¬ 
quired, and /icl can instead be calculated by comparing 
the number of SNe la hosted by cluster galaxies to the 
number that appear hostless and belong to the IC stellar 
population. 

This was first done with the Wise Observatory Optical 
Transients Search (WOOTS) by Gal-Yam et al. (2003), 
who found that 2 of their 7 SNe la in rich, low-redshift 
galaxy clusters appeared to be hostless, which implied 
that /icl ~ 20%. A similar measurement was made with 
the Sloan Digital Sky Survey (SDSS) by McGee & Balogh 
(2010), who found that 19 of the 59 SNe la in low-redshift 
galaxy groups appeared to be hostless, which implied 
that /icl ~ 50% for smaller-scale structures (galaxy 
groups have a total mass ~ 10% that of rich clusters). 
Most recently, the Multi-Epoch Nearby Cluster Survey 
(MENeaCS) from the Canada France Hawaii Telescope 
(CFHT) identified four apparently hostless SNe la in low- 
redshift massive galaxy clusters (Sand et al. 2011), which 
implied that f IC L = O.ie/g^g. 

This usefulness of SNe la as tracers of the ICL depends 
on their being truly hostless. Stacked imaging from the 
past surveys left 0.03 — 0.3% (WOOTS), 3% (SDSS), and 
0.05 — 0.1% (MENeaCS) of the stellar mass in cluster 
galaxies below the detection thresholds. These values are 
dependent on the logarithmic faint-end slope of the lu¬ 
minosity function, ad (he. the Schechter function). The 


above results are based on ad = —1.0, which is true for 
field galaxies (Blanton et al. 2003), but ad may be higher 
in rich galaxy clusters (e.g., Milne et al. 2007). Adopt¬ 
ing ad = —1.5 as an upper limit, the MENeaCS survey 
estimated that < 2% of the stellar mass in dwarf cluster 
galaxies below detection thresholds. If the SNIa occur¬ 
rence rate per unit mass is equivalent in all cluster stel¬ 
lar populations - detected galaxies (~ 82%), undetected 
dwarf galaxies (< 2%), and intracluster stars (~ 16%) - 
then for MENeaCS we expect to find that < 2% of all the 
discovered cluster SNela (23) are hosted by faint dwarf 
galaxies (~ 0.5 SNela). In this work we use deep HST 
imaging at the locations of the 4 MENeaCS IC SNela 
to further lower the fraction of undetected mass in faint 
galaxies to 0.2% (0.003-0.007% with ad = —1-0), analyze 
previously undetected objects in the vicinity of 2 SNela, 
and discuss how our results affect the ability of IC SNe la 
to measure /icl at higher redshifts. 

1.2. SNIa Delay Times 

One path towards understanding the progenitor sce¬ 
nario of SNe la is to constrain the range of possible ages 
via the SNIa delay time distribution (DTD; the time be¬ 
tween star formation and explosion). While Type II and 
Ib/c supernovae - the core collapse of massive stars - are 
associated with young stellar populations, SNela occur 
with an explosion rate proportional to both galaxy stel¬ 
lar mass and star formation rate (e.g., Mannucci et al. 
2005; Scannapieco & Bildsten 2005; Sullivan et al. 2006). 
This indicates that SNela are associated with both old 
and young stellar populations. The current best mea¬ 
surements of the SNIa DTD indicate that some SNIa 
progenitors are quite old, >2 Gyr (e.g. Maoz, Mannucci 
& Nelemans 2014), which implies long-lived progenitors 
and/or that the timescales for mass transfer or merger 
are long. 

The conclusion that at least some SNIa require pro¬ 
genitors that are > 2 Gyr old relies on the SNIa rate 
in rich cluster galaxies, where the majority of the stellar 
population is old (e.g. Sand et al. 2012). However, there 
is evidence that low levels of star formation are present 
in elliptical galaxies (e.g., Yi et al. 2005; Suh et al. 2011; 
Graham et al. 2012). Could it be that all SNIa are ac¬ 
tually from younger stellar populations? Probably not - 
for example, Graham et al. (2012) show that the SNIa 
DTD result is robust to this small amount of young stars. 
Even so, direct confirmation of a SNIa progenitor in a 
stellar population of purely old stars would strengthen 
and support the late-time DTD constraints on the pro¬ 
genitor scenario. 

A suitable environment for this test is the population of 
intracluster (IC) stars in rich galaxy clusters: the colors 
and luminosities of IC red giant stars in Virgo show that 
the IC stellar population is comprised of stars > 2 Gyr 
old (Durrell et al. 2002), and theoretical models also 
suggest the ICL is comprised of old stars (e.g. Sommer- 
Larsen et al. 2005; Purcell et al. 2007). All of the 
apparently hostless cluster SNe discovered to date are of 
Type la; a lack of intracluster core collapse supernovae 
lends credence to the idea that the ICL is composed of 
an older stellar population. In this work, we use deep 
HST imaging to show that at least two of the apparently 
hostless SNela discovered by MENeaCS truly belong to 
the IC stellar population of old stars. 


1.3. SNIa Rates in Faint Hosts 

Over the past ten years, an increasing number of SN 
surveys have employed an unbiased, wide-held search 
strategy instead of targeting massive galaxies. This has 
lead to the discovery that some types of luminous SNe 
have a higher explosion rate in dwarf hosts (e.g. Neill et 
al. 2011; Lunnan et al. 2014), attributed to high star 
formation rates providing more progenitor stars and/or 
lower metallicity leading to more luminous SNe, or per¬ 
haps an elevated rate of binary star formation. As most 
SNela occur in massive galaxies, it is difficult to assess 
whether they might also have an enhanced rate in dwarf 
galaxies because of the large statistical uncertainty on 
the rate due to the relatively low number of SNela in 
low-mass hosts (e.g. Quimby et al. 2012). SN surveys of 
rich galaxy clusters are an efficient way to search many 
low-mass galaxies at once, both due to a higher sky den¬ 
sity of galaxies and the (putative) upturn in the faint-end 
slope of the cluster luminosity function (e.g., Milne et al. 
2007; Yamanoi et al. 2007). In this work we find that 
one IC SNIa may be hosted by a dwarf cluster galaxy, 
and discuss the implication of this for SN la rates in faint 
hosts. 

GCs are another potential very faint host for IC SNe la. 
No SN has ever been confirmed to be hosted by a GC, 
although they are theoretically predicted to have a SN la 
occurrence rate 1 — 10 x that of elliptical galaxies due to 
dynamical interactions that lead to more white dwarfs in 
binary systems in GCs (e.g., Ivanova et al. 2006; Sliara & 
Hurley 2006; Pfahl et al. 2009). Non-detections of GCs 
at the sites of ~ 45 low-redshift SNela in archival HST 
images have constrained the potential rate enhancement 
to < 42 x the rate in elliptical galaxies (Voss & Nelemans 
2012; Washabaugh & Bregman 2013). In this work we 
find that one IC SN la may be hosted by a GC, and dis¬ 
cuss the implication of this for SNIa rates in GCs. Since 
GCs are also comprised mainly of old (5-13) Gyr stellar 
populations, confirming a GC-hosted SNela would also 
meet our science goal of confirming SNe la with long (> 2 
Gyr) delay times. 

1.4. Paper Overview 

In this work, we use the Hubble Space Telescope (HST) 
Advanced Camera for Surveys (ACS) to obtain deep im¬ 
ages in filters F606W and F814W at the locations of 4 
apparently hostless SNela from MENeaCS (Sand et al. 
2011). These are the deepest images, and the largest 
single survey sample, of IC SNela locations in rich clus¬ 
ters ever obtained and analyzed in this way. In Section 
[2] we present our original CFHT and new HST observa¬ 
tions and discuss our image processing and photometric 
calibrations. In Section [3] we analyze our deep stacks of 
HST ACS imaging: we characterize the faint sources in 
the vicinity of the SNela, derive our point-source lim¬ 
iting magnitudes, and determine the amount of cluster 
stellar luminosity remaining below our detection thresh¬ 
olds. In Section [4] we discuss the implications of these 
results with respect to SNIa progenitor ages, the rates 
of SNela in faint cluster objects, and the future use of 
SNela as tracers of the ICL. We conclude in Section [5j 
All dates are given in UT and a standard flat cosmology 
of Dm = 0.3, Da = 0.7 is assumed throughout. 

2. OBSERVATIONS 


Here we describe our past observations with the Mega- 
Cam imager (Boulade et al. 2003) at the Canada-France- 
Hawaii Telescope (CFHT), and our new deep imaging 
taken 4-5 years later with the Hubble Space Telescope 
HST Advanced Camera for Surveys (ACS; Holland et al. 
1998). 

2.1. CFHT MENeaCS 

The Multi-Epoch Nearby Cluster Survey (MENeaCS) 
monitored 57 low-redshift (0.05 < z < 0.15) massive 
galaxy clusters from 2008-2010 with the MegaCam in¬ 
strument at the Canada-France-Hawaii Telescope. A to¬ 
tal of 23 cluster SNela (Sand et al. 2012) and 7 cluster 
SNell (Graham et al. 2012) were discovered. The sur¬ 
vey strategy, spectroscopic follow-up, detection efficien¬ 
cies, and the derivation of SN rates from MENeaCS are 
presented in Sand et al. (2012), and the four intracluster 
SNe la discovered by MENeaCS are presented in Sand et 
al. (2011). For the 4 IC SNela, in Table |T) we list their 
MENeaCS identifier, coordinates, spectroscopic redshift, 
and type as determined by the Supernova Identification 
(SNID) software package (Blondin & Tonry 2007). As 
described in Sand et al. (2011), these SNela are all 
within 1 Mpc of the brightest cluster galaxy and within 
3000 km s -1 of the cluster redshift, thereby confirming 
they occurred in the cluster. All four were more than 
5 effective radii (i.e., > 5x the half-light radius) from 
the nearest potential host galaxy in the CFHT images. 
Deep stacks were made from SN-free CFHT survey im¬ 
ages, and implanted simulated point sources were used 
to constrain the limiting magnitude of any possible host 
galaxy. These limits are listed for each IC SN la in Ta¬ 
ble [I] along with an upper limit on the fraction of light 
in low-mass cluster galaxies below our detection limit ( 
/(< L m i n )). Assuming a faint-end slope for the cluster 
luminosity function of ad = —1.5, the fraction of stellar 
mass below the detection thresholds of our CFHT deep 
stacks is <2% (Sand et al. 2011). 

2.2. HST ACS Imaging 

To assess whether these four apparently hostless SNe la 
truly belonged to the intracluster stellar population we 
used the HST ACS to obtain deep images at their loca¬ 
tions. We waited until >3 years after peak brightness 
to obtain these images to avoid contamination from the 
SN itself (this is discussed in Section ET5l) . To either rule 
out or classify a faint object at the SN position we re¬ 
quired an imager that offers both sensitivity and high 
resolution. Although the Wide Field Camera 3 (WFC3) 
has slightly better resolution, we chose the ACS for its 
higher throughput and increased sensitivity beyond 4000 
A because most objects in galaxy clusters are red. At 
our target position on the CCD we selected aperture 
“WFC” because amplifier “B” on WFC1 has the low¬ 
est read noise. We restricted our telescope orientations 
to avoid possible stellar diffraction spikes coming near 
our target coordinates. Because a photometric color is 
necessary to classify a detected object we used two wide 
filters: F606W (H-band) and F814W (/-band). We di¬ 
vided our integration times into multiple exposures to 
remove cosmic rays and used small dithers to mitigate 
the effect of hot pixels, and were able to obtain all ob¬ 
servations of a given cluster in a single visit. We provide 
a summary of our HST observations in Table [2] 


To make the deepest possible stacks of our /7ST+ACS 
data, we use the Astrodrizzle software to median-combine 
FLC images provided by that STScI pipeline (Gonzaga 
et al. 2012). FLC files come fully reduced, drizzled, 
and corrected for charge transfer efficiency by the STScI 
pipeline. We use the WCS astrometry during image com¬ 
bination. We do this for the F606W and F814W images 
separately to make a deep stack for each filter that is free 
of cosmic rays, and then also create a single sum-combine 
image for the deepest possible stack. 

2.3. HST ACS Photometric Calibrations 

To obtain the apparent magnitude of a source, inf , 
where / represents filter F606W or F814W, we start 
with the raw magnitude, m/ >raw , for which we use 
MAG_AUTO from Source Extractor (Bertin & Arnouts 
1996). We add a small PSF correction for point sources, 
derived from the simulated sources used for the limiting 
magnitudes (~ —0.1, see Section HTTl) . This is not usually 
necessary with MAG_AUTO but we find it is required for 
the non-Gaussian PSF of the HST ACS (the Tiny Tim 
PSF; Krist et al. 2011). We also apply the zeropoint, 
zj = —2.5 log Pf — 21.1, where Pf is the PHOTFLAM 
header keyword, representing the flux of a source with 
constant Fa which produces a count rate of 1 electron 
per second (Pf has units of erg cm -2 s -1 A -1 ). For 
our images, ZF 606 W = 26.66 and 3 f 814W = 26.78 mag. 
To transform from the natural system and obtain the 
apparent magnitude, mp, where F is the Johnson V or 
Cousins I filters, we use this equation from Sirianni et al. 
(2005; see their Equation 12): 

TUf = TUf raw + Co f + (Ci f x TCOL) + (C2 F x TOOL 2 ) 

(!) 

where TCOL is the color in the targeted system; in our 
case TCOL = my — mi- We use the synthetic coefficient 
values from Sirianni et al. (2005; see their Table 22), as 
listed in our Table [3] 

3. ANALYSIS 

In this section we analyze our deep HST ACS images 
at the locations of our 4 intracluster SNela from ME¬ 
NeaCS. We determine our limiting point-source magni¬ 
tudes in Section m and derive the fraction of cluster 
luminosity remaining below our detection thresholds in 
Section 13.21 We present and characterize any objects 
found near the SN la coordinates in each cluster in Sec- 
tion l3.31 quantify the likelihood that detect ed ob jects are 
the result of a chance alignment in Section 13.41 and rule 
out the possibility of observing the evolved companion 
star or lingering SN la emission in Section 13.51 

3.1. Point Source Limiting Magnitudes 

In order to determine the limiting magnitude of our 
images, we plant 5000 fake stars in each of our stacked 
F606W and F814W images. These fake stars have ap¬ 
parent magnitudes 26.0 < m < 30.0, with more stars 
at fainter magnitudes, and an appropriate PSF from the 
Tiny Tim model PSF generator (Krist et al. 2011). We 
ensure that simulated stars are only planted in regions of 
low surface brightness in order to mimic the locations of 
our intracluster SNe la. We run Source Extractor on the 
images with the simulated population of point-sources; 


TABLE 1 

CFHT Data for MENeaCS Intracluster Supernovae' 


MENeaCS 

Identifier 

SN Coordinates 

RA, Dec 

BCG Offset 
(kpc) 

Redshift 

(SNID) 

Spectral 

Type 

CFHT Detection Limit 
Mg, M r 

f L-min ) 
a d = -1.5 

Abelll650_9_13_0 

12:59:01.33, -01:45:51.68 

468 

0.0836 

la-norm 

-12.47, 

-13.04 

0.0172 

Abell2495_5_13_0 

22:50:26.33, +10:54:41.70 

148 

0.0796 

la-norm 

-11.72, 

-12.37 

0.0127 

Abell399_3_14_0 

02:57:26.41, +12:58:07.63 

616 

0.0613 

la-norm 

-12.54, 

-12.56 

0.0138 

Abell85 6 08 0 

00:42:02.39, -09:26:58.00 

595 

0.0617 

Ia-91bg 

-11.15, 

-11.68 

0.0091 


a From Table 1 in Sand et al. (2011). 


TABLE 2 

HST+ACS Imaging Data for MENeaCS Intracluster Supernovae 


MENeaCS 

Identifier 

Orbits 

Observation 
Date (UT) 

Exposure Time (s) 
F606W F814W 

HST Detection Limits (mag) 
F606W F814W Mr 

f Lmin) 
OLd = -1.5 

Abelll650_9_13_0 

5 

2013-01-25 

3870 

6120 

28.58 

28.88 

-9.7 

0.0029 

Abell2495_5_13_0 

4 

2013-10-04 

2984 

5970 

28.64 

28.90 

-9.8 

0.0026 

Abell399_3_14_0 

2 

2013-12-04 

1800 

1980 

28.30 

28.45 

-10.1 

0.0035 

Abell85_6_08_0 

1 

2013-09-18 

1800 


28.40 


-9.2 

0.0020 


1 

2014-09-05 


1800 


28.49 




TABLE 3 

Photometric Transformations 


Natural 

Filter 

Target 

Filter 

Co 

C i 

c 2 

TCOL 

my — mi 

F606W 

Johnson V 

26.394 

0.153 

0.096 

< 0.4 

F606W 

Johnson V 

26.331 

0.340 

-0.038 

> 0.4 

F814W 

Cousins I 

25.489 

0.041 

-0.093 

< 0.1 

F814W 

Cousins I 

25.496 

-0.014 

0.015 

> 0.1 


the relevant detection parameters are given in Table [|] 
We visually verify that these parameters are returning 
all and only real sources in the images. In Figure [T| we 
plot our detection efficiency (i.e. the fraction of objects 
recovered) as a function of apparent magnitude for the 
F606W and F814W images of each cluster. The “limit¬ 
ing magnitude” is defined as the magnitude at which the 
detection efficiency drops to 50%, and is given for each 
cluster and filter in the plot legend of Figure [TJ 

In Section 13.31 below, we visually identify objects be¬ 
low the official limiting point source magnitude at the 
locations of our SNela in Abell 2495 and Abell 399. We 
find that these objects are only detected by Source Ex¬ 
tractor if the threshold is lowered to lcr, which also de¬ 
tects many peaks in the background noise and produces a 
source catalog with large uncertainties in their apparent 
magnitude. This is why our official limiting magnitude 
- which needs to be robust because we use it to deter¬ 
mine the fraction of cluster stellar luminosity below our 
detection thresholds in Section 13.21 - is slightly brighter 
than some of the sources discussed in Section 13.31 The 
caveat here is that brighter, but more extended, objects 
may fall below our detection threshold also - but most of 
the faint cluster objects will be point-like (dwarf galaxies 
and GCs). 

3.2. Fraction of Undetected Cluster Light 

To determine the fraction of cluster light remaining 
below our point-source limiting magnitudes, we follow a 
similar method to that presented by Sand et al. (2011) 
for our CFHT deep-stack images, the result of which is 
reproduced in the last column of Table [l] 


TABLE 4 

Source Extractor Parameters 
for the Detection Efficiency 


Parameter 

Value 

DETECT.MINAREA 

3 

DETECT.THRESH 

2 

SEEING.FWHM 

0.08 

BACK.SIZE 

256 

BACK.FILTERSIZE 

5 

BACK.TYPE 

AUTO 

BACKPHOTO.TYPE 

LOCAL 


Point Source Limiting Magnitudes 



Simulated Apparent Magnitude 


Fig. 1.— Detection efficiencies for simulated p oint s ources in our 
HST ACS images, as described in the text (§ 13.111 . Filter and 
cluster are represented by color and symbol respectively, as shown 
in the plot legend. The magnitude at which our detection efficiency 
falls to 50% is considered our limiting magnitude. 


The absolute I?-band luminosity function for the 
nearby Virgo cluster is modeled by Trentham & Tully 
(2002) with two components, a Gaussian for the bright 

















end and a Schechter function (see their Equation 2): 

N(M) = N g e- {M - M o )2/(2 ^ ] 

+ N d (10 [ - 0 ^ M - Md)] ) ad+1 e ~ lol ~°' i(M ~ Mdn , ( 2 ) 

where N (M) is the number density of galaxies per square 
Mpc per magnitude bin, M is the /?-band absolute 
magnitude, and N g = 17.6, M g = —19.5, a g = 1.6, 
Nd = 3N g , M d = —18.0, and a d = —1.03. However, the 
faint-end slope is known to steepen with redshift (e.g., 
Khochfar et al. 2007), and values down to a d ~ —1.5 
have been measured for the Coma cluster (Milne et al. 
2007) and most recently for Abell 85 ( a d « —1.6; Agulli 
et al. 2014). In order to make a robust upper limit on the 
amount of cluster light below our detection efficiencies, 
we use a d = —1.5 from here on. 

In Sand et al. (2012) the total r'-band luminosity is 
calculated for all MENeaCS clusters. We convert this 
to i?-band using the conversion factors from Blanton & 
Roweis (2007), R = r - 0.0576 - 0.3718((r- i) - 0.2589). 
Because most of the cluster light is from old stellar popu¬ 
lations, we use the typical SDSS color of elliptical galax¬ 
ies, r — i ~ 0.4, from Eisenstein et al. (2001). We inte¬ 
grate the galaxy luminosity function down to the CFHT 
point-source limiting magnitudes listed in Table Q] from 
Sand et al. (2011), and then normalize to the total R- 
band luminosity for each cluster. 

We convert the absolute i?-band luminosity function 
for each cluster into apparent V- and /-band using the 
cluster’s redshift, the elliptical galaxy template spec¬ 
trum from Kinney et al. (1996), and the line-of-siglrt 
Galactic extinction for each cluster: Ay,Ai 650 = 0.047, 
Av,A2495 = 0.211, Ay,A399 = 0.467, Ay,A85 = 0.103, 
4li,Ai650 = 0.026, A.j ; A 2495 = 0.116, Ai : A 399 = 0.256, and 
Ai a 85 = 0.05 7 (Schlafly & Finkbeiner 2011; Schlegel, 
Finkbeiner & Davis 1998). We convert the V and /-band 
luminosity function into the HST ACS natural system fil¬ 
ters using the transformations of Sirianni et al. (2005), 
as described in Section 12.31 In Table [2j we report the 
point-source detection limit in absolute Mr magnitudes, 
and the fraction of the cluster’s /?-band stellar luminosity 
remaining below this limit. For all clusters we find that 
< 0.2% of the stellar mass in cluster galaxies remains 
below our official point-source limiting magnitude. This 
result is discussed further in Section |4] 

3.3. Potential Hosts in the Deep HST-ACS Images 

We use the IRAF tasks GEOMAP and GEOTRAN to 
co-register our CFHT images to the new, deeper HST- 
ACS images. In Figures [2] [5] we show the results, side by 
side, for comparison. The CFHT images are comprised of 
two 120 second r' exposures, and contain the SN la - we 
do not use our SN-free deep stacks here, because we need 
the SN’s coordinates in the co-registered frames. The 
HST images are our deepest stacks, the sum-combined 
F606W and F814W filtered images. Green circles mark 
the position of the SN in each image, with a radius equal 
to 3 x the positional uncertainty of the SN la. These posi¬ 
tional uncertainties, listed in Table 0 are a combination 
of Source Extractor’s uncertainty in the PSF centroid 
for the SN in the co-registered CFHT image (using win¬ 
dowed output parameters), and the error in GEOMAP’s 
transformation between images. The dashed cyan lines 


TABLE 5 

Supernova Positional Uncertainty® 


SN la 

PSF Center 13 

GEOMAP c 

Total d 

Error 

Cluster 

Ax 

Ay 

Ax 

Ay 

Ax 

Ay 

Radius 6 

Abell 

(pix) 

(pix) 

(pix) 

(pix) 

(pix) 

(pix) 

H 

1650 

4.57 

4.57 

0.31 

0.13 

7.52 

7.62 

0.23" 

2495 

6.55 

6.31 

1.81 

2.00 

6.42 

6.46 

0.34" 

399 

6.15 

6.00 

0.63 

1.56 

7.76 

7.82 

0.31" 

85 

7.39 

7.21 

2.57 

3.52 

10.41 

10.15 

0.40" 


a All pixels are HST ACS pixels (0.05" pixel x ). 
b Source Extractor’s uncertainty on the SN’s coordinates in 
the CFHT image post-transformation with GEOTRAN (i.e., the 
square root of the variance, the second moment of the barycenter). 

c GEOMAP’s uncertainty in mapping the CFHT image to the 
HST image (i.e., a systematic, the output x r -ms and y r ms values). 
d Added in quadrature. 
e Average in x and y , x0.05 // pixel -1 . 

enclose nearby objects (sizes chosen to guide the eye). 
Image information such as the cluster name, UT date of 
acquisition, filters, scale bar, and compass are shown in 
yellow along the bottom. These images all have a 30" x 
30" field of view. 

As in Sand et al. (2011), we use the dimensionless pa¬ 
rameter R to identify whether nearby objects could be 
considered as potential hosts of the SNela. This param¬ 
eter is defined in the Source Extractor manual as: 

R — C xx X r Gyyy r T C X yX r y r (3) 

where C xx , C yyi and C xy are object ellipse parameters, 
x r = XSN-Xgai and y r = ysN~y ga l, and R ~ 3 describes 
the isophotal limit of the galaxy (see also Sullivan et al. 
2006). A supernova is typically only classified as “host¬ 
less” if R > 5, but depending on the surface brightness 
profile of the galaxy, a significant amount of the stellar 
mass may reside beyond this radius (e.g., > 10% for the 
potential host of the SN in Abell 399, as determined in 
Sand et al. 2011). 

3.3.1. Abell 1650 

In Figure [2j we see that the location of the SN la is 
truly devoid of objects. We ran Source Extractor with 
very relaxed parameters and still recovered no sources in 
this area. Of the three nearby objects enclosed by dashed 
circles in Figured the SN location is R > 15 away. Our 
intracluster SNIa in Abell 1650 therefore appears to be 
truly hostless. 

3.3.2. Abell 2495 

In Figure [3] we see several sources near the location 
of the SNIa that were not apparent in the CFHT deep 
stacks presented in Sand et al. (2011), but none are 
within the positional uncertainty of the SNIa. We ran 
Source Extractor with very relaxed parameters and still 
recovered no sources within the positional uncertainty. 
In Figure [6] we show this region in detail, and identify 
nearby sources A, B, C, and D. Sources A and B are 
unlikely to be physically associated, as the SNIa is R ~ 
17 away from them. We discuss objects C and D in turn. 

Object C is likely a part of object D, which is clumpy 
and extended, as shown by the contour plot in the right- 
hand image of Figure [6] However, object C is identified 
by Source Extractor as an independent source at the lcr 






Fig. 2.— Co-registered CFHT (left) and HST (right) images for the SNIa in Abell 1650, as described in the text of Section 13.31 



Fig. 3. — Co-registered CFHT (left) and HST (right) images for the SNIa in Abell 2495, as described in the text of Section 13.31 


level in the F814W image, with tofsitw = 29.0±0.2 mag. 
In the F606W image it is not officially detected by Source 
Extractor, but it is just visible to the eye, and with aper¬ 
ture photometry we estimate it to be 771F606W = 29.8±0.2 
mag. In both filters, object C falls below our 50% detec¬ 
tion efficiency for Abell 2495 (see Figure[T]in Section 15711) . 
The SN’s location is R > 5 away from object C, which 
argues against a physical association. Object C is redder 
than the red sequence of Abell 2495, as shown in Figure 
0 and so it is unlikely to be a cluster dwarf galaxy. The 
scenarios in which object C is lingering emission from the 


SN la or a chance alignment are discussed in Sections 13.41 
and 13.51 Ultimately we find it unlikely that object C is 
the host galaxy. 

Object D is an extended source with semi-major and 
semi-minor axes of A = 0.22" and B = 0.053", which 
leads to SN offsets of R = 4.9 and 4.1, in F606W and 
F814W respectively. As R < 5, this SNIa cannot be 
considered truly hostless unless we can show that object 
D is unlikely to be a cluster member. In the natural 
system of HST ACS, object D is TOF 606 W = 25.56 and 
WF 814 W = 25.75 mag, and has a color = —0.20 mag. 




Fig. 4.— Co-registered CFHT (left) and HST (right) images for the SNIa in Abell 399, as described in the text of Section 13.31 



Fig. 5. — Co-registered CFHT (left) and HST (right) images for the SNIa in Abell 85, as described in the text of Section 13.31 


As shown in Figure [7] it consistent with a red sequence 
cluster galaxy. We use Sirianni et al. (2005) to convert 
this photometry into Landolt filter system and find that 
object D is my ~ 25.6 and m; « 24.5 mag. We apply 
the distance modulus of Abell 2495 (/j, « 37.8 mag) and 
find that intrinsically, object D is My « —12.2 and Mi ~ 

— 13.3 mag. This is brighter than the limiting magnitudes 
quoted for the CFHT deep stack of Abell 2495, but those 
limits are for point-like sources and object D is extended 

- in fact, it is just barely visible as an extended source 
in Figure 3 of Sand et al. (2011). 


Object D is clumpy and has an ellipticity of 0.7, which 
is higher than the ellipticity of the bright red sequence 
galaxies we identify in Abell 2495. Morphologically, ob¬ 
ject D resembles an inclined disk galaxy - can we use 
the disk scale length to assess whether it may be a blue 
spiral galaxy at higher redshift? Disk galaxies are gen¬ 
erally well fit by an exponential function for the flux in¬ 
tensity as a function of radius, I(R) = Ioe~ R ^ Rd , where 
Rd is the characteristic disk scale length. We estimate 
Rd < 0.4", with no attempt to de-project or account for 
inclination. If object D is a cluster member, this corre- 








Fig. 6.— The HST +ACS image from the deepest stack of both filters F606W and F814W for the intracluster SNIa in Abell 2495. The 
green circle marks the position of the SN la, with a radius equal to 1 x the positional uncertainty from Table [5] In cyan we identify nearby 
objects A, B, C, and D (we cannot verify whether C is a part of D), discussed in the text. For objects A and B, the dashed cyan circle has 
r = 0.24", representing 3x the PSF FWHM of ACS (0.08"). For the dashed ellipse around object D we use the parameters determined by 
Source Extractor for the F814W image: 3x the semi-major and -minor axes, A = 4.4 and B = 1.1 pixels (A = 0.22" and B = 0.053"), and 
a position angle of 6 = 74 degrees (in the x-y plane). In the image at right, we include flux contours in magenta to highlight the connection 
between objects C and D and the lack of low-significance sources within the SN’s positional error circle. 
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Fig. 7. — Color-magnitude diagram in the natural HST ACS filter 
system for Abell 2495. Black points represent all sources extracted 
with an apparent magnitude brighter than the 50% dete ction ef¬ 
ficiencies in both F606W and F814W, as discussed in § EH and 
shown in Figure IT1 Black points surrounded by red indicate the 
galaxies used in the linear fit to the red sequence (red lines). The 
blue square and circle represent objects C and D respectively (see 
FigurelGl). 

sponds to Rd ~ 0.6 kpc, which is roughly appropriate for 
a disk galaxy. In Section [374] we also discuss the relatively 
large probability that this is a chance alignment, but for¬ 
mally we cannot exclude the possibility that object D is 
a cluster member and the host of the SN in Abell 2495. 

3.3.3. Abell 399 

In Figure Q] we see a large elliptical galaxy near the lo¬ 
cation of the SNIa. Sand et al. (2011) describes how the 
SN is > 7 R from the large elliptical, and how the SN’s 
and galaxy’s line-of-siglit velocities differ by 3<r, indicat¬ 
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ing they are not associated. We use these HST images to 
re-evaluate the parameter R from Equation [3j which de¬ 
scribes the SN’s offset in terms of the galaxy’s isophotes. 
We find that i? ~ 6 in both F606W and F814W, re¬ 
establishing the SN as independent of the stellar popu¬ 
lation of this large neighbor galaxy. 

In Figure [4] we see a small, faint source at the center 
of the green circle marking the SN la’s position, labeled 
object F in Figure [5] Object F is detected with low sig¬ 
nificance, but with coincidence in both the F606W and 
F814W images. It is a point-like source, with a FWHM 
of ~ 3 pixels, and is offset by < 1 pixel from the SNIa 
location. Object F is therefore very likely to be physi¬ 
cally associated with the SN la in Abell 399 - the small 
probability of a chance alignment with a foreground star 
or background host is discussed in Section EH In Fig¬ 
ure [5] we also label the next-nearest object E; it is a 
clumpy extended source, but we find it far less likely to 
be physically associated with the SN and so we focus our 
attention on object F. 

Object F is TOF 606 W = 28.7 ± 0.3 and tof 814 W = 
28.7 ± 0.2 mag, below our limiting magnitudes as shown 
in Figure [T) Its photometry is consistent with a red se¬ 
quence galaxy in Abell 399, which is very similar to the 
red sequence shown for Abell 2495 in Figure [3 We use 
Sirianni et al. (2005) to convert this photometry into the 
Landolt filter system and find that object F is my ~ 28.8 
and to/ « 27.4 mag. We apply the distance modulus of 
Abell 399 (/jl ss 37.2 mag) and find that intrinsically, ob¬ 
ject F is My ~ —8.4 and M/ Ri —9.8 mag, and has a 
color of V — I ~ 1.4. This is also well matched to the 
expected magnitude and color of a bright GC: the lu¬ 
minosity function for GCs is a Gaussian that peaks at 
Mi rs —7.4 and has er Ri 1.2, and GC colors span from 
My — Mi r i 0.7 1.5 (e.g., West et al. 2011). Although 
GCs have been found distributed between galaxies in rich 










Fig. 8.— The HST -\-ACS image from the deepest stack of both 
filters F606W and F814W for the candidate intracluster SN la in 
Abell 399; a zoom-in of Figure [4] 


clusters (e.g. Peng et al. 2011; West et al. 2011), their 
density drops with clustercentric radius and we would 
not expect a significant number of intracluster GCs at 
the location of this SNIa, 616 kpc from the BCG. How¬ 
ever, the presence of the nearby elliptical galaxy (Figure 
a makes the GC hypothesis more likely. 

In Section S3 we show that the SNIa rate enhance¬ 
ments in dwarf galaxies or GC implied by the nature of 
object F are within all existing theoretical and observa¬ 
tional limits. If we assume that the rate per unit mass 
of SNela is not much larger in dwarf galaxies vs. GC 
(or vice versa) - which would put a strong prior on the 
nature of object F - we can estimate the expected sur¬ 
face densities of dwarf red sequence galaxies and red GCs 
at this location in order to assess which is more proba¬ 
ble. Based on our extrapolation of the cluster luminosity 
function for Abell 399 (Section [S3) , we estimate there to 
be ~ 1.5 xlO 9 L 0 in dwarf cluster galaxies at Mr > —10 
mag. That is about 100-500 dwarf galaxies of 5-10 x 10 6 
L 0 within a clustercentric radius of ~ 1 Mpc, or 1-5 
xlO -4 dwarfs kpc -2 . This would apply at any location 
in the cluster, and so applies for the location of object F. 
To assess whether satellites of the nearby elliptical galaxy 
might raise the predicted surface density at this location, 
we use the radial distribution of low-redshift, low-mass 
satellite galaxies presented by Prescott et al. (2011; see 
their Figure 5). The location of the SN in Abell 399 is 
R !=s 10 kpc from the center of the nearby elliptical; at 
this radius Prescott et al. (2011) find that the surface 
density of satellites is ~ 7 x 10 -6 kpc -2 . The caveat here 
is that they consider only isolated primary galaxies, and 
so their results represent an upper limit on the radial dis¬ 
tribution we could expect in rich clusters. However, as 
this is much lower than what we expect from the cluster 
luminosity function we conclude that object F is unlikely 
to be a satellite galaxy of the nearby elliptical. 

We can estimate the surface density of GC from the 


nearby elliptical galaxy (Figure 3), which has absolute 
magnitudes of My ss —19.5 and Mi ss —20.7, and a stel¬ 
lar mass of ~ 7.6 x 10 9 M 0 (Bell et al. 2003). A galaxy of 
this stellar mass has between 10-100 GC within R < 50 
kpc, and the surface density radial distribution for GC 
is N(r) = N 0 r~ 2A GC kpc -2 (Zaritsky et al. 2015); 
~ 75% of the GC are internal to 10 kpc (the distance 
of the SN in Abell 399 from the nearby elliptical). Nor¬ 
malizing the radial distribution to a total of 10-100 GC 
gives N Ki 0.3 —3 x 10 -2 GC kpc -2 . Based on the GC lu¬ 
minosity function of West et al. (2011) ~ 84% of all GC 
are below our limiting magnitudes, and so the observable 
surface density at 10 kpc is lowered to 0.5 — 5 x 10 -3 
GC kpc -2 . This is higher than our estimated surface 
density for cluster dwarf galaxies, suggesting it is more 
likely to see a GC at the location of object F than a 
cluster dwarf. However, as our lower limit for the GC 
surface density is equal to our upper limit on the dwarf 
galaxy surface density the estimates are not significantly 
different enough to make a robust claim to the nature of 
object F. The caveat here is that we have used a radial 
distribution for all GC, but the population of red GCs 
has a significantly shorter radial extent than blue GCs - 
in fact, for cluster ellipticals the radial distribution of red 
GCs appears to be truncated near the effective radius of 
the parent galaxy (e.g., Brodie & Strader 2006). 

In summary, we find that the SNIa in Abell 399 is 
likely physically associated with object F. Without any 
constraints from expected SN la rates, we find that object 
F is less likely to be a dwarf galaxy than a GC from the 
nearby elliptical. We discuss the probability of chance 
alignments in Section 3U and the possibility that object 
F is lingering emission from the SN la in Section 13.51 We 
discuss the implications of object F for SNIa rates in 
dwarf galaxies and GCs, and whether these implications 
provide a prior on the nature of object F, in Section lj.ll 
and the impact of object F on using SNIa as tracers of 
the ICL at high-redshift in Section ITT?! 

3.3.4. Abell 85 

In Figure [5] we see that the location of the SNIa ap¬ 
pears devoid of objects. For the nearest source, labeled 
object G in Figure 3 the SN is R ~ 13 away. We 
ran Source Extractor with very relaxed parameters, but 
found that the few sources detected are consistent with 
noise peaks (i.e., they had very low significance, their de¬ 
tections in F814W only and were not coincident with any 
F606W low-significance detections, and were not visually 
confirmed in the F606W+F814W stack). We conclude 
that the SN la in Abell 85 is truly hostless. 

3.4. Probabilities of Random Line-of-Sight Alignments 

In our discussion of the nature of object F in Abell 
399 we estimated the surface densities of faint cluster 
objects such as dwarf red sequence galaxies and GCs to 
be ~ 5 x 10 -4 objects per kpc 2 . The probability of a faint 
cluster object appearing randomly within the positional 
uncertainty of our IC SNela is negligible, < 0.0003%. 
The presence of object F is not a chance alignment with 
a cluster object. 

We used TRILEGAL (Girardi et al. 2005) to simulate 
a foreground star population in the directions of our four 
fields to a limiting magnitude of my = 29. The prob¬ 
ability of a star appearing randomly within 0.2 " of our 



IC SN la coordinates is <0.02%. We judge that it is ex¬ 
tremely unlikely that object F at the location of the SN 
in Abell 399 is a foreground red star. 

We used the Hubble Deep Field catalogs (Williams et 
al. 1996) to simulate a population of faint objects be¬ 
tween 25.0 < WF 814 W < 29.5. This is the magnitude 
range in which we detect sources in our HST images but 
not our CFHT images. The new sources detected near 
the SN la locations in Abell 2495 and Abell 399 fall in this 
magnitude range. We find the probability of a faint field 
object randomly being within 0.31" and 0.34", the po¬ 
sitional uncertainties on the locations of the SNela, are 
~ 2.0% and ~ 2.5%, respectively. In our original work 
~ 2% of the cluster stellar mass in faint dwarf galaxies 
was below our detection limit (Sand et al. 2011). Assum¬ 
ing the SN la occurrence rate per unit mass is equivalent 
in all populations (high- and low-mass galaxies, and in¬ 
tracluster stars), this means ~ 2% of the all MENeaCS 
SNe la should be hosted by undetected, faint dwarf galax¬ 
ies (~ 0.5 out of 23 SNela). We now see that this is 
approximately the same chance as finding a background 
galaxy within the positional uncertainty. However, ob¬ 
ject F in Abell 399 is not just within the positional un¬ 
certainty, but appears within 0.05" (1 ACS pixel) of the 
SN location. The probability of random alignment with 
a background object within 0.10" is just 0.3%. On the 
other hand, if we increase the radius to 1.0" (i.e. the dis¬ 
tance encompassing nearby objects for Abell 2495), the 
probability of chance alignment increases to ~ 18%. 

As a final, alternative estimate we use our own source 
catalogs for all four HST- ACS fields and find that the 
fraction of our imaged area within R < 5 of an object is 
~ 3%. While it remains very unlikely that object F is a 
chance alignment with the SNIa in Abell 399 (Figure [5]), 
we cannot say the same for the objects near the SNIa 
location in Abell 2495 (Figure [Gj. 

3.5. Limits on the SNela and/or its Companion 

A normal SNIa has faded ~ 4 magnitudes by ~ 100 
days after peak brightness. After this time the decline 
is set by the decay rate of Co 56 , and the SNIa continues 
to fade at ~ 1 magnitude per 100 days in /-band, and 
~ 1.3 — 1.5 mag per 100 days in BVR (e.g. as seen 
for normal twin SNela 201 lfe and 201 lby; Graham et 
al. 2015a). For a normal SNIa such as SN2011fe, the 
intrinsic brightness at ~ 1000 days is V ~ — 5 magnitudes 
(e.g., Kerzendorf et al. 2014; Graham et al. 2015b). 
Most of the flux comes from blue emission features of 
[Fe II] at < 6000 A (e.g., Taubenberger et al. 2015; 
Graham et al. 2015b), and the late-time V — I color of a 
normal SNIa is expected to be ~ 0. After ~ 1000 days, 
the predicted rate of decline for normal SNela is even 
slower (Seitenzahl et al. 2009), and so V ~ —5 mag is a 
conservative upper limit on SNIa brightness after 1000 
days. 

Over-luminous SNela that resemble SN 1991T have 
been observed to decline more slowly. For example, 
SN 1991T itself was V ~ —10 at ~ 600 days (Cappel- 
laro et al. 1997) and SN 2000cx had a F-band slope 
of ~ 0.65 mag per 100 days at ~ 700 days after peak 
brightness (Sollerman et al. 2004), so a SN 1991T-like 
event could be F < — 7 at 1000 days. However, none 
of the MENeaCS SNela were spectroscopically similar 
to SN 1991T, and are furthermore unlikely to belong to 


this subclass because 91T-like SNela are associated with 
younger stellar populations (Howell et al. 2009). Al¬ 
though less is known about the late-time decline of sub- 
luminous SNela, SN 1991bg itself was already F « — 6 
by ~ 600 days after peak brightness (Turatto et al. 
1996). Our intracluster SNIa in Abell 85 was classified as 
SN 1991bg-like, but no object is detected at its location 
in the HST ACS imaging. 

Could we detect the emission from a shocked compan¬ 
ion star? For normal SNela, theoretical predictions for 
a non-degenerate companion shocked by the SN la ejecta 
include an increase in temperature and luminosity, up to 
10 3 to 10 4 L 0 by 1 10 years after explosion (Pan et al. 
2013; Shappee et al. 2013). Such a companion would be 
blue, and have V > —4 mag, which is well below our lim¬ 
iting magnitudes. A light echo is also likely to be blue, 
similar to the color of a SNIa at peak light, and even 
fainter. In order to formally rule out the possibility that 
the objects identified near the SN locations in Abell 2495 
and Abell 399 are lingering emission from the SNIa or 
its binary companion, we discuss each in turn. 

Abell 2495 — In Figure [6] we identify object C as a 
possible point source near the location of the SNIa in 
Abell 2495. The spectrum of this SNIa was obtained 
on 2009-06-18.58 UT at Gemini Observatory as part of 
the MENeaCS follow-up campaign, and we classified it 
as a normal SNIa at ~ +3 months after peak bright¬ 
ness. At the time of our HST ACS images obtained 
on 2013-10-04, the SNIa would be +1569 days old (4.3 
years). At this time, object C is TOF 606 W = 29.8 + 0.2 and 
TOF 814 W = 29.0 ± 0.2 mag, or my ~ 30.0 and mi « 27.8 
mag, in our HST ACS imaging. We apply the distance 
modulus of Abell 2495 {ji « 37.8 mag) and find that 
intrinsically, object C is My « —7.8 and Mi « —10.0 
magnitudes. This is both significantly brighter and red¬ 
der than a normal SNIa is predicted to be at ~ 4 years. 
We conclude that object C is unlikely to be the SNIa or 
an evolved companion star. 

Abell 399 — In Figure [8] we identify object F as a 
point source at the location of the SNIa in Abell 399. 
The classification spectrum of this SNIa, obtained on 
2008-11-28.49 UT at Gemini Observatory as part of the 
MENeaCS follow-up campaign, showed it to be a normal 
SNIa at ~ +2 weeks after peak brightness. In our HST 
ACS images obtained on 2013-12-04, the SNIa would be 
+1832 days old (5 years). This object is My « —8.4 
and Mi —9.8 mag, has a color of F — / s=s 1.4, both 
significantly brighter and redder than a normal SNIa is 
predicted to be at extremely late times. We conclude 
that object F is unlikely to be the SNIa or an evolved 
companion star. 

4. DISCUSSION 

Our analysis of the HST ACS images at the locations 
of our 4 intracluster MENeaCS SNela has shown that 
one is hosted by either a dwarf red sequence galaxy or 
red GC (Abell 399); one is potentially associated with a 
nearby spiral or irregular galaxy consistent with the red 
sequence but also has a relatively high probability of be¬ 
ing a chance alignment (Abell 2495); and two appear to 
be truly hostless (Abell 1650 and Abell 85). We discuss 
the implications of our results for the rates of SNela in 
faint cluster hosts in Section |4~T1 and for the use of SNela 
as tracers of the ICL in Section 14.21 


4.1. Implications for SNIa Rates in Clusters 

In Section [3] we found that the SN la in Abell 399 was 
likely hosted by the faint object F, and that this source 
is consistent with being either a cluster dwarf galaxy or 
a GC. Here we consider the implications of both scenar¬ 
ios on the rate of SNela in dwarf galaxies and GCs, and 
whether established SNIa rates (or limits) in these pop¬ 
ulations can constrain the physical nature of object F. 

4.1.1. Dwarf Galaxies 

If object F in Abell 399 is a dwarf galaxy, does this 
imply a significantly enhanced SN la rate in faint cluster 
galaxies? The rate per unit mass in a population, 1Z, 
is expressed by 1Z = C x N/M, where C is a detection 
efficiency, N is the number of SNe la, and M is the mass 
in the population. As described in Sand et al. (2011), 
our original CFHT deep stacks left < 2% of the mass 
in faint cluster galaxies undetected, but we now believe 
that population has hosted 1 SN la. We can estimate the 
implied relative rate per unit mass in the faintest 2% of 
cluster galaxies with the following equation: 

7 ^2% _ C 2% x N 2 %/M 2 % , 4 . 

1Z C x N/M ' 1 ’ 

Sand et al. (2011) describes how a small difference 
in MENeaCS detection efficiencies between hosted and 
hostless SNela are introduced by two effects: (1) it is 
more difficult to detect transients on top of a host galaxy 
(even with difference imaging techniques), and (2) the 
spectroscopic follow-up coverage for the hostless popu¬ 
lation was slightly more extensive than for the hosted 
SNe (they were run under separate proposals). Together, 
this difference works out to be C 2 % = 1.2 C. Sand et 
al. (2012) present that the number of SNela hosted by 
all cluster galaxies within 1 Mpc is N = 11, and so 
with N 2 % = 1 we find that 1Z 2 %/1Z « 5.5. Repeating 
this calculation using only red sequence cluster galaxies 
yields a similar rate enhancement because the number of 
SNe la in red sequence members within 1 Mpc is Nr s = 6 
(i.e., 0.5 AT), the stellar mass in red sequence galaxies is 
Mrs ~ 0.5M, and C RS = C. 

Ultimately this potential rate enhancement by a fac¬ 
tor of ~ 5 is quite uncertain, as it is based on just one 
SN la and an indirect estimate of the amount of mass in 
faint cluster galaxies. As introduced in § 11.11 and 13.41 
by assuming the SNIa rate per stellar mass is equal in 
all cluster populations we estimated that the expectation 
value for the number of MENeaCS SNela in dwarf hosts 
is ~ 0.5. For a more restrictive estimate of the expecta¬ 
tion value, we limit to red sequence galaxies within the 
clustercentric radius of 1 Mpc used above. Assuming the 
rate in bright cluster red sequence galaxies (~ 50% of 
the stellar mass hosting TVrs = 6 SNela) is the same as 
that in dwarf red sequence galaxies (~ 2% of the stellar 
mass), the expectation value is ~ 0.24 SNela in dwarf 
hosts. With Poisson statistics the probability of observ¬ 
ing > 1 SNIa given this expectation value is 0.16, and 
the la uncertainties our estimated rate enhancement are 
5.5(_ 2 4 7 5 (Gelirels 1986) - consistent with no rate enhance¬ 
ment. Furthermore, if the luminosity function is steeper 
than expected there could be more than 2% of the stellar 
mass residing in such faint red galaxies. We plan to use 
our deep HST images to constrain the faint-end slope of 


the cluster luminosity function in later work, but con¬ 
sider it beyond the scope of this paper. 

The potential SNIa rate enhancement by a factor of 
~ 5 in the faintest dwarf galaxies is not too large to be 
unphysical, as the SNIa rate is known to vary by factors 
of that size (and more) with host galaxy properties such 
as the specific star formation rate (e.g. Mannucci et al. 
2005; Scannapieco & Bildsten 2005; Sullivan et al. 2006; 
Smith et al. 2012). This potential rate enhancement is 
also not so large that we expect to have already observed 
it in wide-held sky surveys such as SDSS. For example, 
Smith et al. (2012) find that the SNIa rate per unit 
mass decreases as a function of host stellar mass, but 
their lowest stellar mass bin is ~ 5 x 10 8 M e and the 
uncertainty on its rate is a factor of ~ 2-3. Furthermore, 
they mention that ~ 2% of the SNela in their sample 
have undetected host galaxies. Current and future wide- 
held surveys such as the Palomar Transient Factory and 
the LSST should be able to improve the SNIa rate in 
faint dwarf galaxies (e.g., Conroy & Bullock 2015). 

4.1.2. Globular Clusters 

If object F in Abell 399 is a GC, it would be the hrst 
confirmed GC to host a SNIa. As GCs are purely old 
stellar populations (> 2 Gyr), such an association would 
also be direct confirmation that SNela progenitors in¬ 
clude truly old star systems. The fraction of galaxy stel¬ 
lar mass in GC systems is ~ 2 x 1CV 3 for elliptical galax¬ 
ies (Harris et al. 2009; Peng et al. 2008; Zaritsky et 
al. 2015). Based on this and the total stellar mass in 
red sequence galaxies in MENeaCS clusters, we estimate 
that > 1O 12 M 0 in GCs has been surveyed by MENeaCS. 
If object F is a GC and has produced a SN la, it implies 
a rate ~25 times higher than the rate in our cluster red 
sequence galaxies from Sand et al. (2012). This is on 
a similar scale to theoretically predicted enhancements 
due to dynamical interactions in the dense stellar en¬ 
vironments of GCs (Pfahl et al. 2009), and below the 
current limits placed by non-detections of GC at the lo¬ 
cations of SNela in archival HST images (< 42x; Voss 
& Nelemans 2012; Washabaugh & Bregman 2013). 

Would such an enhanced rate in GCs have already been 
noticed? Potentially not. It would imply that < 5% of 
the SNela in ellipticals ar e host ed by their GCs, but 
as we discussed in Section 13.3.31 the radial distribution 
of GCs follows the galaxy light profile and drops off at 
~ 51?. It is entirely conceivable that GC-hosted SNe have 
simply been assigned as belonging to the parent galaxy. 
This might lead to a relative rate enhancement in the 
outer regions of ellipticals, which could be measured and 
attributed to GCs. However, such an effect would be 
difficult to confidently measure for two reasons. First, 
there is a detection bias of SNe being easier to find when 
they are not embedded deep in the host galaxy. Second, 
a stellar population originating in GCs - able to produce 
an enhanced rate of SN la - may have previously released 
into the halos (or bulges) of galaxies from GCs due to col¬ 
lisions with clouds or other GCs, winds from supernovae 
in the GC, tidal forces for GC on elliptical orbits about 
their host, evaporation of stars during interal GC relax¬ 
ation, and/or dynamical friction (Fall & Rees, 1985). It 
is conceivable that these issues combined could conspire 
to blur the signal of a SN la rate enhancement in GCs in 
the radial distribution of SNela. 






4.1.3. Rates Summary 

We find that the implied SNIa rate enhancements in 
either dwarf galaxies or GCs do not conflict with exist¬ 
ing observations or theoretical predictions, but also do 
not provide a means to constrain the nature of object F. 
Additionally, if either of the 2 SNela in Abell 1650 and 
Abell 85 - which appear hostless in our deep HST images 
- were in fact hosted by the ~ 0.2% of the stellar mass 
in small galaxies that remains below our limiting magni¬ 
tudes, then by Equation [4] this indicates a rate enhance¬ 
ment by a factor of ~ 55 in the faintest dwarf galaxies. 
This is conspicuously high and would have been noticed 
in wide-field surveys. We therefore conclude that these 2 
SNela were truly hosted by the population of intraclus¬ 
ter stars stripped from their host galaxy and residing in 
the cluster’s gravitational potential. 

4.2. Implications for IC SNela as Tracers of the ICL 

Sand et al. (2011) reported that of the 23 cluster 
SNela discovered by MENeaCS, 4 had no apparent host 
galaxy in deep CFHT images that left just ~ 2% of the 
total cluster stellar mass in undetected faint galaxies. 
With our deep HST imaging we find that 1 out of the 
23 cluster SNe la, ~ 4%, is hosted by a faint point-source 
(object F in Abell 399). Finding ~ 4% of the SNela 
hosted by ~ 2% of the stellar mass is not surprising, but 
here we take a closer look at how /icl is derived in order 
to confirm the utility of SNe la as tracers of the ICL at 
higher redshift. 

The fraction of intracluster light, /icl, is calculated 
by dividing the number of hostless SNela by the total 
number of cluster SNela (hosted+hostless) discovered 
by MENeaCS (Sand et al. 2011). To do this, the de¬ 
tection efficiencies must be equal for hosted and hostless 
SNela; in other words, a survey must apply the same 
discovery and spectroscopic classification constraints to 
both populations, or be able to account for any bias in 
the pipeline. This was true for all MENeaCS SNela ex¬ 
cept the IC SN in Abell 2495, which was preferentially 
observed with Gemini spectroscopy despite being fainter 
than the magnitude limit applied to follow-up of ME¬ 
NeaCS SNe. For this reason, the hostless SN in Abell 
2495 was not included in the calculation of the fraction 
of IC stellar mass in Sand et al. (2011). In order to 
combine only the same physical regions of each cluster, 
they limit to a radius < i? 2 oo (he., the virial radius); all 
four apparently hostless SNela are within i? 2 oo- Using 
the remaining 3 apparently hostless SNela as intraclus¬ 
ter, and the 13 hosted SNela within i? 2 oo, Sand et al. 
(2011) measure /icl = 0.16±o;J|. 

Given the proximity of the SNIa in Abell 399 to a 
large nearby galaxy (see Section [3.3.311 . they also repeat 
this calculation assuming that this SN la was hosted, and 
find /icl = 0.11^0 or hi light of the faint object F at 
the location of the SNIa in Abell 399, we can now say 
the latter is the more accurate measurement of /icl • Al¬ 
though this distinction may not seem significant because 
the difference between these two /icl measurements is 
within their relatively large statistical errors, future facil¬ 
ities such as the Large Synoptic Survey Telescope (LSST) 
will generate bigger sample sizes and have smaller un¬ 
certainties and a better understanding of the fraction of 
apparently hostless SNela will be needed in this regime. 


LSST itself will be able to provide this because its deep 
stack images have a projected detection limit of r ~ 27.5 
mag0. 

Until then, we can only caution that future surveys 
will have to consider that 25-30% of apparently hostless 
SNela might not belong to the IC stellar population. Is 
that too large for hostless SNe la to be scientifically use¬ 
ful tracers of /icl at higher redshifts? Some numerical 
simulations indicate that /icl grows with cosmological 
time, and by z ~ 0 is ~ 2x larger than at 2 ~ 0.4 (e.g. 
Murante et al. 2007); others find that most IC stars 
are stripped at z > 1, in which case /icl would remain 
constant since then (Puchwein et al. 2010). Direct mea¬ 
surements of this low surface brightness component have 
shown that /icl 25% at z < 0.1 and ~ 10% at z ~ 0.2 
(e.g. Gonzalez et al. 2007; Zibetti et al. 2005), but HST 
imaging of0.4<2<0.8 clusters has found no evolution 
hi /icl since z < 0.8 (Guennou et al. 2012). We there¬ 
fore surmise that assessments of /icl from SNela will 
require uncertainties of < 30% in order to compare with 
some theoretical models and the direct surface bright¬ 
ness measurements. This sounds discouraging, but there 
is hope: below, we suggest that the best way to improve 
this uncertainty and use hostless SNela as high-redshift 
ICL tracers is to constrain the SNIa occurrence rate in 
faint hosts. 

Dwarf galaxies represent a small fraction of the total 
cluster stellar mass, which is constrained by measure¬ 
ments of the galaxy luminosity function. If their SNIa 
occurrence rate is equal to that in elliptical galaxies, then 
measurements of /icl can account for the contamination 
from apparently hostless SNe la in dwarf hosts. However, 
in the preceding section we show that if object F is a clus¬ 
ter dwarf, the SN la occurrence rate in the faintest ~ 2% 
of cluster galaxies could be up to ~ 5x higher than in 
elliptical galaxies. If confirmed, the number of SNela in 
faint dwarfs could be up to half of all apparently hostless 
SNela. 

Compared to dwarf galaxies, GC represent an even 
smaller fraction (< 0.002) of the total stellar mass in 
clusters - much less than the fraction of intracluster stars 
(~ 0.16). If their SNIa occurrence rate is equal to that 
in elliptical galaxies, accidentally including the very few 
SNe la in GC at large radial offsets from their host as part 
of the ICL will have a negligible effect on measurements 
of ficL- (Most SNela in GC will be associated with el¬ 
liptical galaxies, because GC have a radial distribution 
similar to that of stars.) On the other hand, if the SNIa 
occurrence rate in GC is 25 x that in elliptical galaxies 
- as implied if object F is a GC - then up to ~ 5% of 
all cluster SNela, and up to ~ 30% of the apparently 
hostless cluster SNela, may actually be associated with 
GCs. 

In this work we have shown that ~ 75% of all ap¬ 
parently hostless SNela are truly intracluster, and that 
contamination is only a significant problem if the SNIa 
occurrence rate in low-mass galaxies or GC is enhanced. 
This issue can be resolved by SNIa rates analyses from 
low-redshift, wide-field surveys such as the Palomar 
Transient Factory, SDSS, or LSST, combined with deep 
imaging for a larger sample of low-redshift apparently 
hostless SNIa. Such an effort should sufficiently con- 

5 http://lsst.org/lsst/science/science_portfolio 




strain the SNIa rate in dwarfs and GCs such that ex¬ 
tremely deep imaging will not be required to confirm all 
apparently hostless cluster SNe la, and facilitate their use 
as tracers of the ICL to higher redshifts. 

5. CONCLUSION 

We have presented deep HST +ACS images at the lo¬ 
cations of 4 IC SNe la in rich galaxy clusters, obtained 
>3 years after explosion. This is the largest single-survey 
sample of IC SNe la in rich clusters, and these data are 
the deepest images ever obtained at the locations of IC 
SNe la, lowering the amount of stellar mass left unde¬ 
tected in cluster galaxies from ~ 2% to just ~ 0.2% 
(< 0.005% if we assume a shallow faint-end slope for 
the galaxy luminosity function, ad = —1.0). We have 
confirmed that the 2 SNe la in Abell 1650 and Abell 85 
are hostless, and truly belong to the intracluster stellar 
population of stars stripped from their host galaxy and 
residing in the cluster potential. This indicates that at 
least some SNIa progenitors have truly old progenitor 
stars (> 2 Gyr). We could not rule out that the SNIa in 
Abell 2495 was hosted by a nearby disk galaxy that has a 
magnitude, color, and size consistent with cluster mem¬ 
bership, but also found a relatively high probability that 
this is a random association. We have shown that the 
SN la in Abell 399 was very likely hosted by a faint red 
point-like source that has a magnitude and color consis¬ 
tent with both dwarf red sequence galaxies and red GCs. 
Our statistical analysis of the expected surface densities 
has shown that a dwarf galaxy is less likely at that loca¬ 
tion than a GC, due to the presence of a nearby elliptical 
galaxy. We have demonstrated that the rate enhance¬ 
ments in dwarfs or GCs implied by this new faint host are 
plausible under current observational constraints, and we 
do not reject either hypothesis. We have also explicitly 
ruled out the possibility that we could observe extremely 
late-time emission from the SNe la themselves or a pos¬ 
sible shocked companion. 

Our discovery of a faint host galaxy for 1 of the 4 


Agulli, I., Aguerri, J. A. L., Sanchez-Janssen, R., Barrena, R., 
Diaferio, A., Serra, A. L. and Mendez-Abreu, J. 2014, MNRAS, 
444, 34 

Bell, E. F., McIntosh, D. H., Katz, N. and Weinberg, M. D. 2003, 
ApJS, 149, 289 

Bertin, E. and Arnouts, S. 1996, A&AS, 117, 393 
Blanton, M. R. et al. 2003, ApJ, 592, 819 
Blanton, M. R. and Roweis, S. 2007, AJ, 133, 734 
Blondin, S. and Tonry, J. L. 2007, ApJ, 666, 1024 
Brodie, J. P. and Strader, J. 2006, ARA&A, 44, 193 
Boulade et al. 2003, in Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference Series, Vol. 4841, Society of 
Photo-Optical Instrumentation Engineers (SPIE) Conference 
Series, ed. M. Iye & A. F. M. Moorwood, 72-81 
Cappellaro, E., Mazzali, P. A., Benetti, S., Danziger, I. J., 

Turatto, M., della Valle, M. and Patat, F. 1997, A&A, 328, 203 
Conroy, C. and Bullock, J. S. 2015, ApJL, [arXiv: 1504.04015] 
DeMaio, T., Gonzalez, A. H., Zabludoff, A., Zaritsky, D. and 
Bradac, M. 2015, MNRAS, 448, 1162 
Durrell, P. R. and Ciardullo, R. and Feldmeier, J. J. and Jacoby, 
G. H. and Sigurdsson, S. 2002, ApJ, 570, 119 
Eisenstein, D. J., Annis, J., Gunn, J. E., Szalay, A. S., Connolly, 
A. J. et al. 2001, AJ, 122, 2267 
Fall, S. M. and Rees, M. J. 1985, ApJ, 298, 18 
Feldmeier, J. J., Ciardullo, R., Jacoby, G. H. and Durrell, P. R. 
2004, ApJ, 615, 196 


SNe la that appeared to be hostless is a potential prob¬ 
lem for measuring the evolution in the fraction of intra¬ 
cluster light across cosmic time, because /icl will need to 
have uncertainties < 25% in order to distinguish between 
models or compare with independent measurements. For 
this reason we argue that hostless SNe la in rich clusters 
should be used to measure /icl with caution, and that 
it would be best if deep imaging for a larger sample of 
low-redshift, apparently hostless SNe la were obtained in 
order to better constrain the rate in dwarf galaxies and 
GCs. 

This work is based on observations made with the 
NASA/ESA Hubble Space Telescope and obtained from 
the Mikulski Archive for Space Telescopes at the Space 
Telescope Science Institute, which is operated by the 
Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS 5-26555. These observa¬ 
tions are associated with program #12937 and partially 
funded by a NASA grant. We thank the HST staff for 
their assistance with this program. 

This work is also based in part on observations ob¬ 
tained with MegaPrime/MegaCam, a joint project of 
CFHT and CEA/DAPNIA, at the CanadaFranceHawaii 
Telescope (CFHT) which is operated by the National Re¬ 
search Council (NRC) of Canada, the Institut National 
des Science de lUnivers of the Centre National de la 
Recherche Scientifique (CNRS) of France, and the Uni¬ 
versity of Hawaii. We remain grateful to the operators of 
the CFHT queue without whom the original MENeaCS 
project would not have been possible. 

MLG’s position in the supernova research group at 
U.C. Berkeley is supported by Gary & Cynthia Bengier. 
MLG and DJS were supported during HST proposal 
preparation by the Las Cumbres Observatory Global 
Telescope Network. CJP is supported by the Natural 
Sciences and Engineering Research Council of Canada. 
We thank our anonymous reviewer for their constructive 
comments and useful suggestions. 


Gal-Yam, A., Maoz, D., Guhathakurta, P. and Filippenko, A. V. 

2003, AJ, 125, 1087 
Gehrels, N. 1986, ApJ, 303, 336 

Girardi, L., Groenewegen, M. A. T., Hatziminaoglou, E. and da 
Costa, L. 2005, A&A, 436, 895 

Gonzaga, S., Hack, W., Fruchter, A., Mack, J., eds. 2012, The 
DrizzlePac Handbook. (Baltimore, STScI) 

Gonzalez, A. H., Zabludoff, A. I. and Zaritsky, D. 2005, ApJ, 618, 
195 

Gonzalez, A. H., Zaritsky, D. and Zabludoff, A. I. 2007, ApJ, 666, 
147 

Graham, M. L. et al. 2012, ApJ, 753, 68 
Graham, M. L. et al. 2015, MNRAS, 446, 2073 
Graham, M. L., Nugent, P. E., Sullivan, M., Filippenko, A. V., 
Cenko, S. B., Silverman, J. M., Clubb, K. I. and Zheng, W. 
2015, MNRAS, [arXiv: 1502.00646] 

Guennou, L. et al. 2012, A&A, 537, 64 

Harris, W. E., Kavelaars, J. J., Hanes, D. A., Pritchet, C. J. and 
Baum, W. A. AJ, 137, 3314 
Holland, F. C. et al. 1998, SPIE, 3356, 234 

Howell, D. A., Sullivan, M., Brown, E. F., Conley, A., Le Borgne, 
D. et al. 2009, ApJ, 691, 661 
Howell, D. A. 2011, Nature Communications, 2, 350. 

Ivanova, N., Heinke, C. O., Rasio, F. A., Taam, R. E., Belczynski, 
K. and Fregeau, J. 2006, MNRAS, 372, 1043 



Kerzendorf, W. E., Taubenberger, S., Seitenzahl, and I. R., 
Ruiter, A. J., 2014, ApJ, 796, 26 
Kinney, A. L., Calzetti, D., Bohlin, R. C., McQuade, K., 
Storc.hi-Bergmann, T., Schmitt, H. R. 1996, ApJ, 467, 38 
Khochfar, S., Silk, J., Windhorst, R. A., and Ryan, R. E., Jr. 
2007, ApJ, 668, 115 

Krist, J. E., Hook, R. N. and Stoehr, F. 2011, SPIE Conference 
Series, 8127, 0 

Lunnan, R., Chornock, R., Berger, E., Laskar, T., Fong, W., Rest, 
A., Sanders, N. E., Challis, P. M., et al. 2014, ApJ, 787, 138. 
Neill, J. D., Shara, M. M. and Oegerle, W. R. 2005, ApJ, 618, 692 
Neill, J. D., Sullivan, M., Gal-Yam, A., Quimby, R., Ofek, E., 
Wyder, T. K., Howell, D. A., Nugent, P. et al. 2011, ApJ, 727, 
15 

Maoz, D., Mannucci, F., and Nelemans, G. 2014, ARA&A, 52, 107 
Mannucci, F., Della Valle, M., Panagia, N., Cappellaro, E., 

Cresci, G., Maiolino, R., Petrosian, A. and Turatto, M. 2005, 
A&A, 433, 807 

McGee, S. L., and Balogh, M. L. 2010, MNRAS, 403, L79 
Milne, M. L., Pritchet, C. J., Poole, G. B., Gwyn, S. D. J., 

Kavelaars, J. J., Harris, W. E. and Hanes, D. A. 2007, AJ, 133, 
177 

Montes, M. and Trujillo, I. 2014, ApJ, 794, 137 
Murante, G., Giovalli, M., Gerhard, O., Arnaboldi, M., Borgani, 
S. and Dolag, K. 2007, MNRAS, 377, 2 
Pan, K.-C., Ricker, P. M., Taam, R. E. 2013, ApJ, 773, 49 
Peng et al. 2008, ApJ, 681. 197 
Peng E. et al. 2011, ApJ, 730, 23 

Pfahl, E., Scannapieco, E. and Bildsten, L. 2009, ApJ, 695, 111 
Prescott, M. et al. MNRAS, 417, 1374 
Puchwein, E., Springel, V., Siiacki, D. and Dolag, K. 2010, 
MNRAS, 406, 936 

Purcell, C. W.. Bullock, J. S. and Zentner, A. R. 2007, ApJ, 666, 
20 


Quimby, R. M., Yuan, F., Akerlof, C., Wheeler, J. C., Warren, 

M. S. 2012, AJ, 144, 177 
Sand, D. J. et al. 2011, ApJ, 729, 142 
Sand, D. J. et al. 2012, ApJ, 746, 163 
Scannapieco, E. and Bildsten, L. 2005, ApJ, 629, 85 
Schlafly E. F., Finkbeiner D. P., 2011, ApJ, 737, 103 
Sc.hlegel D. J., Finkbeiner D. P. and Davis M., 1998, ApJ, 500, 
525 

Seitenzahl, I. R., Taubenberger, S. and Sim, S. A. 2009, MNRAS, 
400, 531 

Shara, M. M. and Hurley, J. R. 2006, ApJ, 646, 464 
Sirianni, M. et al., PASP, 117, 1049 

Shappee, B. J., Kochanek, C. S., Stanek, K Z. 2013, ApJ, 765, 150 
Smith, M. et al. 2012, ApJ, 755, 61 
Sollerman, J. (et al. 2004, A&A, 428, 555 
Sommer-Larsen, J., Romeo, A. D. and Portinari, L. 2005, 
MNRAS, 357, 478 

Suh, H., Yoon, S.-c., Jeong, H. and Yi, S. K. 2011, ApJ, 730, 110 
Sullivan, M. et al. 2006, ApJ, 648, 868 
Taubenberger, S. et al. 2015, MNRAS, 448, 48 
Trentham, N. and Tully, R. B. 2002, MNRAS, 335, 3 
Turatto, M., Benetti, S., Cappellaro, E., Danziger, I. J., Della 
Valle, M., Gouiffes, C., Mazzali, P. A. and Patat, F. 1996, 
MNRAS, 283, 1 

Voss, R. and Nelemans, G. 2012, A&A, 539, 77 
Washabaugh, P. C. and Bregman, J. N. 2013, ApJ, 762, 1 
West, M. J., Jordan, A., Blakeslee, J. P., Cote, P., Gregg, M. D., 
Takamiya, M. and Marzke, R. O. 2011, A&A, 528, 115 
Williams, R. E. et al. 1996, AJ, 112, 1335 
Yamanoi, H. et al. 2007, AJ, 134, 56 
Zaritsky, D. et al. 2015, ApJ, 799, 159 

Zibetti, S., White, S. D. M., Schneider, D. P. and Brinkmann, J. 
2005, MNRAS, 358, 949 



